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By means of a new mod(N) invariant operator basis, s-parametrized phase-space functions associ- 
ated with bounded operators in a finite-dimensional Hilbert space are introduced in the context of 
the extended Cahill-Glauber formalism, and their properties are discussed in details. The discrete 
Glauber-Sudarshan, Wigner, and Husimi functions emerge from this formalism as specific cases of 
s-parametrized phase-space functions where, in particular, a hierarchical process among them is 
promptly established. In addition, a phase-space description of quantum tomography and quantum 
teleportation is presented and new results are obtained. 



I. INTRODUCTION 

The first proposal of a unified formalism for quasiprobability distribution functions in continuous phase space has 
its origin in the seminal works produced by Cahill and Glauber y| . Since then a huge number of papers have appeared 
in the literature covering a wide range of practical applications in different physical systems modeled by means of 
infinite-dimensional Hilbert spaces Q. In particular, the phase-space description of some important effects in 
quantum mechanics, such as interference, entanglement, and dccohcrence, has opened up astounding possibilities for 
the comprehension of intriguing aspects of the microscopic world Q. However, if physical systems with a finite- 
dimensional space of states are considered, then the quasiprobability distribution functions are described by a set of 
discrete variables defined over a finite lattice [[IJLLjLLj] In this sense, Opatrny et al were the 

first researchers to propose a unified approach to the problem of discrete quasiprobability distribution functions in the 
literature. Basically, they used a discrete displacement-operator expansion to introduce s-parametrized phase-space 
functions associated with operators defined over a finite-dimensional Hilbert space. Furthermore, the authors showed 
that the discrete Glauber-Sudarshan, Wigner, and Husimi functions are particular cases of s-parametrized phase- 
space functions and depend on the arbitrary reference state whose characteristic function cannot have zero values. 
It is worth mentioning that the dependence on the right choice of the reference state and the associated problems 
with the mod(iV) invariance of the discrete displacement operators represent two important restrictions inherent to 
their approach which deserve to be carefully investigated. Nowadays, beyond these fundamental features, discrete 
quasiprobability distribution functions in finite-dimensional pha se spaces have potential applications for quantum- 
state tomography 0, , quantum teleportation 0, 0, 120I |2 ll . phase-space representation of quantum computers 
[22| , open quantum systems |2^| , quantum information theory |24| , and quantum computation [23 ■ 

The main aim of this paper is to present a consistent formalism for the quasiprobability distribution functions 
defined over a discrete ./V 2 -dimensional phase space, which is based upon the mathematical fundamentals developed 
in |26| . First, we review important topics and introduce new properties concerning the mod(TV) invariant opera- 
tor basis which leads us not only to define a parametrized phase-space function in terms of the discrete s-ordered 
characteristic function, but also to discuss some characteristics inherent to the extended Cahill-Glauber formalism 
for finite-dimensional spaces. The restriction on the right choice of the reference state is overcome in this approach 
through the vacuum state established by Galetti and de Toledo Piza 8] , whose analytical properties were extensively 
explored in |lOj |. Consequently, the discrete Glauber-Sudarshan, Wigner, and Husimi functions are well-defined in 
the present context and represent specific cases of s-parametrized phase-space functions describing density operators 
associated with physical systems whose space of states is finite. In addition, we also establish a hierarchical order 
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among them through a smoothing process characterized by a discrete phase-space function that closely resembles the 
role of a Gaussian function in the continuous phase-space. In this point, it is worth emphasizing that our ab initio 
construction inherently embodies the discrete analogues of the desired properties of the Cahill-Glauber approach. 
Next, we apply such discrete extension into the context of quantum information processing, quantum tomography, 
and quantum teleportation in order to obtain a phase-space description of some topics related to unitary depolarizers, 
discrete Radon transforms, and generalized Bell states. In particular, we attain new results within which some of 
them deserve to be mentioned: (i) we show that the symmetrized Schwinger operator basis introduced in Q can be 
considered a unitary depolarizer; (ii) we establish a link between measurable quantities and s-ordered characteristic 
functions by means of discrete Radon transforms, which can be used to construct any quasiprobability distribution 
functions defined over a iV 2 -dimensional phase space; and finally, (iii) we present a quantum teleportation protocol 
that leads us to reach a generalized phase-space description of the physical process discussed by Bennett et al • 

This paper is organized as follows. In section II we present some basic properties inherent to the new discrete 
mapping kernel which allow us to define a parametrized phase-space function in terms of a discrete s-ordcred charac- 
teristic function. Following, in section III we show that the extended Cahill-Glauber formalism not only introduces 
new mathematical tools for the analysis of finite quantum systems, but also can be applied in the context of quan- 
tum information processing, quantum tomography, and quantum teleportation. Moreover, we also employ a slightly 
modified version of the scattering circuit to measure any discrete Wigner function in the phase-space representation. 
Finally, section IV contains our summary and conclusions. 



II. THE MAPPING KERNEL 



There is a huge variety of probability distribution functions defined in continuous quantum phase-spaces whose range 
of practical applications in physics covers different areas and scenarios 0, 0] . For example, the well-known Cahill- 
Glauber formalism 1] provides a general mapping technique of bounded operators which permits, in particular, to 
define a generalized probability distribution function F^ s \q,p) = Tr[T^(g,p)p] associated with an arbitrary physical 
system described by the density operator p. In this approach, the mapping kernel (hereafter h = 1) 

TW(?.P) = / ^^exp[i(^-p' g )]DW( g ',p') (1) 

is defined as a Fourier transform of the parametrized operator 

D !, V,P') = exp[( S /4)( 9 ' 2 + p' 2 )]B(q',p') (2) 

where T>(q',p') — exp[i(p'Q — q'P)] is the usual displacement operator written in terms of the coordinate and mo- 
mentum operators satisfying the Weyl-Heisenberg commutation relation [Q,P] = il, and s is a complex parameter. 
Thus, for s = — 1,0, +1 the generalized probability distribution function leads to the so-called Husimi, Wigner and 
Glauber-Sudarshan functions, respectively. Besides, these functions present specific properties and correspond to 
different ordered power-series expansions in the annihilation and creation operators of the density operator: the 
Husimi function !K(q,p) is infinitely differentiable and it is associated with the normally ordered form; the Wigner 
function W(q,p) is a continuous and uniformly bounded function, it can take negative values and corresponds to 
the symmetrically ordered form; and finally, the Glauber-Sudarshan function 7 > (q,p) is highly singular; it does not 
exist as a regular function for pure states and it corresponds to the antinormally ordered form. After this condensed 
review of the Cahill-Glauber formalism for the quasiprobability distribution functions, we will establish the discrete 
representatives of these functions in an iV 2 -dimensional phase space. 



A. The new mod(iV) invariant operator basis 

Let us introduce the symmetrized version of the unitary operator basis proposed by Schwinger |27j| as 

S(r ? ,0 = ^expfe^UW« (3) 

where the labels r\ and £ are associated with the dual coordinate and momentum variables of a discrete iV 2 -dimensional 
phase space. Consequently, these labels assume integer values in the symmetrical interval [—£,£], with £ — (N — l)/2. 
A comprehensive and useful compilation of results and properties of the unitary operators U and V can be found 
in reference ^(j > since the initial focus of our attention is the essential features exhibited by 10 ■ Note that the set 
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of iV 2 operators {S(?7, £)}ri,£=-£,...,t constitutes a complete orthonormal operator basis which allows us, in principle, 
to construct all possible dynamical quantities belonging to the system 27]. Thus, the decomposition of any linear 
operator O in this basis is written as 

l 

o= Yu &(v,OS(v,0 (4) 

with the coefficients 0(?7,£) given by Tr[St(?7, £)0]. It must be stressed that this decomposition is unique since the 
relations S^(ry,^) = S(— 77, — £) and Tr[S^(rj, £)S(r/, £')] = 5^5^ are promptly verified. The superscript [N] on the 
Kronecker delta denotes that this function is different from zero when its labels are mod(iV) congruent. 
The new mod(N) invariant operator basis recently proposed in pfij . 



27ri 

i7r*(»7,6JV)- — (w + £v) 



S (s) (r?,0. (5) 



n,t=-e 

is defined by means of a discrete Fourier transform of the extended mapping kernel 

sW(r ? ,0 = [3C(r?,e)]- s S(r ? ,0 

where the extra term 3C(r],t;) can be expressed as a sum of products of Jacobi theta functions evaluated at integer 
arguments |28j. 

%( m £) = {2 [tf 3 (a|ia)i? 3 (0|4ia) + i? 4 (0|ia)z9 2 (0|4ia)]} _1 {i} 3 (Trari\ia)tf 3 (na£\ia) + t? 3 (7ra77|ia)79 4 (7ra£|ia) exp(i7r?7) 

+ ??4(7ra?7|ia)^3(7ra^|ia) exp(i7r£) + ^4(7ra7y|ia)'i?4(7ra^|ia) exp[i7r(?7 + £ + N)]} (6) 

with a = (2A r ) _1 . As mentioned in [2(|, 3C(t7, £) is a bell-shaped function in the discrete variables (77, £) and equals to 
one for 7/ = £ = 0; in addition, the complex parameter s obeys \s\ < 1. The phase $(r],!;;N) = Nl 1^ — rjl^ — £lf 

is responsible for the mod (AT) invariance of the operator basis if = W/N] being the integral part of a with 
respect to N. This definition stands for the discrete version of the continuous mapping kernel and represents the 
cornerstone of the present approach. 

By analogy with decomposition ipjl. the expansion 

= 1 E 0(- s )( i i,,)TW( ;J ,,) (7) 

H t v=—l 

can also be verified for any linear operator. Here, the coefficients 0^ s ^ (fx, v) = Tr[T^ s ^ (/i, v)0] correspond to a 
one-to-one mapping between operators and functions belonging to an A^ 2 -dimensional phase space characterized by 
the discrete labels /i and v. In particular, if one considers s = — I and O = p in equation lf?Jl. we obtain the diagonal 
representation 



I 1 

- £ nw)\w){»A (8) 



AT 



where 1/) = Tr [T^(fi, v)p\ is the discrete version of the Glauber-Sudarshan function for finite Hilbert spaces, 
and T^O, v) is the projector of discrete coherent-states [l0| . For s = 0, we verify that 



1 ' 

- E W(Ai,i/)G(Ai,i/) (9) 



" AT 

recovers the well-established results in [8j, W(/x, f) = Tr[Gt(/i, i/)p] being the discrete Wigner function and G(/j,, v) the 
mod(A^) invariant operator basis whose mathematical properties were studied in |ldj| . Furthermore, we note that the 
Husimi function in the discrete coherent state representation, 9~C(li, v) — Tr[T^ _1 ^(^t, v)p\, can be promptly obtained 
from equations © or @ by means of a trace operation. Next, we will discuss some properties inherent to the set of 
N 2 operators {T^'(/i, v)}u u =— l with emphasis on establishing a hierarchical process among the quasiprobability 
distribution functions in finite-dimensional spaces. 
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B. Basic properties 

The discrete mapping kernel T^- 1 (fj,, v) presents some inherent mathematical features that lead us to derive a set 
of properties which characterize its algebraic structure. For instance, it is straightforward to show that the equalities 



a) i x: tW( M ,,)=i 

H,v=-l 

(ii) Tt[t« {h,v)\ =1 
(hi) Tr[T(- s W)T( s Vy) 



are promptly verified where, in particular, the third property has been reached with the help of the auxiliary relation 



Tr 



1 

N 



M)+e(^-^)]M3C(^0] 



-(*+«) 



Note that for s = —1, the first pro pert y coincides with the completeness relation of the discrete coherent states (the 
proof of this relation was given in [jij]); the second property simply states that T^ s ^(/Lt, v) has a unit trace. Finally, 
the third property is the counterpart to the orthogonality rule established for the operators S(?y,£). Furthermore, we 
also verify the condition T^ s >(fi, v) — [T( s )(/i, u)]\ which implies that for real values of the parameter s, the discrete 
mapping kernel is Hcrmitian; consequently, the mappings of Hermitian operators in the iV 2 -dimensional phase space 
lead us to obtain real functions. Now, let us establish a hierarchical process among the discrete Glauber-Sudarshan, 
Wigner and Husimi functions. 

The connection between the discrete Glauber-Sudarshan and Wigner functions is reached with the help of equation 
||HJ| through a smoothing process of CP(/i, v), i.e., 



N 



E E(//- M y-^)ovy) 



(10) 



p,' ,v'=-l 



where E(// — /i, v' — v ) = Tr[T(°)(/Lt,^)T(- 1 )(/i', v')\ is expressed by means of a discrete Fourier transform of the 
function 3C(?y,£) - note that E(/i' — fi, v' — v) can be interpreted as a Wigner function evaluated for the discrete 
coherent states labeled by p! and v' . Similarly, the link between discrete Wigner and Husimi functions can also be 
established through equation © as follows: 



1 1 



(ii) 



Therefore, equations and exhibit a sequential smoothing which characterizes a hierarchical process among 
the quasiprobability distribution functions in finite-dimensional spaces, 3 (/i, v) — > — > !K(fi,v). It is worth 
mentioning that 



1 e 



(12) 



establishes an additional relation which allows us to connect both the discrete Husimi and Glauber-Sudarshan func- 
tions without the intermediate process given by "W(n,v), being | (/x, v\(jf, v') | 2 = \0C(/i' — fi, v 1 — v)\ 2 the overlap 
probability for discrete coherent states. Opatrny et al |2( have used a similar formalism in order to establish a set 
of parametrized discrete phase-space functions for finite-dimensional Hilbcrt spaces, where some mathematical pro- 
cedures were introduced to circumvent the condition of mod(iV) invariance of the discrete displacement operators. In 
that approach, the discrete s-parametrized functions basically depend on the arbitrary reference state whose charac- 
teristic function cannot have zero values. Here, we have established a suitable mathematical procedure that allows 
us to overcome some intrinsic problems encountered in [ij] , being the vacuum state defined in [3, H3 as our reference 
state. 
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Next, we present two important properties associated with the trace of the product of two bounded operators and 
the matrix elements (m|T^ s ^(/x, v)\n) in the finite number basis {|n)} n= o ! ... ! jv-i- The first one corresponds to the 
overlap 



i 1 

(iv) Tr(AB) = — ]T A^M^^v) 



where, in particular, for s = 0, the trace of the product of two density operators coincides with the overlap of the 
discrete Wigner functions of each density operator, 



In addition, the mean value of any bounded operator can also be obtained from this property, 



i i 

0}^Tr(Op) = - Q^'W^'W) 



(13) 



being the parametrized function F( s >(jjl, v) defined as the expectation value of the discrete mapping kernel i.e. 



F«(/x,*/) = Tr TW(/i,i/)p 



is) I 



J2 ex p 



ri£=- 



2tti 



sW(t? s o 



(14) 



while S^(?7,£) = Ti^S^^, £)p] represents the discrete s-ordered characteristic function [[J. Note that <I>(?7,£;iV) 
can be discarded in equation (|14fl since the discrete labels rj and £ are confined into the closed interval [—£,£]. In 
fact, this phase will be important only in the mapping of the product of M quantum operators lQj. Besides, for 
s = — 1,0, +1 the parametrized function is directly related to the discrete Husimi, Wigner and Glauber-Sudarshan 
functions, respectively. Hence, the characteristic function can now be promptly calculated for each situation through 
the inverse discrete Fourier transform of the generalized probability distribution function F^(fi, u). 
The second one refers to the nondiagonal matrix elements in the finite number basis 



1 ' 

(v) HT«0i,i/)|n> = - Y, ex P 



2-ki 



[X(r,^)}- s r mn ( m 



with 



r m „(?7,C) = exp ( —fiVZ) Y exp ( ~W (7T l)$w : fc-t,T. 
written in terms of the coefficients Hi 



(15) 



f—l\n 00 
3tt,n = N n ^= 2J eX P 



/3=-oo 



+ ^rP K ) H » [\l ^ 



N 



N 



where N„ is the normalization constant, and H„(z) is a Hermite polynomial. It is easy to show that r mn (?7, £) satisfies 

the relations r m „(0, 0) = Sm} n and Too(r),£) = ^iv^Oi which ar e associated with the orthogonality rule for the 
finite number states and the diagonal matrix element (0|T^ (/x, f)|0) for the vacuum state. Moreover, adopting the 
mathematical procedure established in 29] for the continuum limit, we obtain 



r~T / / i ■ / \ n—m 



(n-m) f lg +iP 



exp 



(?i > m) 



with L^ m) (z) being the associated Laguerre polynomial. Consequently, the nondiagonal matrix elements for \s\ < 1 
take the analytical form 



|TW(g,p)|» 



1 — s V 



(-1-:)" 




7i — m 

t (n—m) 






1 - s 




1 - s 2 



exp 



q 2 +p 2s 
' l-s , 
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This result coincides exactly with that obtained by Cahill and Glauber pj for the mapping kernel (JTJ, since T^ s \/j,, v) 
goes to T( s \q,p) in the limit N — > oo. Following, we will discuss some applications for the generalized probability 
distribution function F^ s \fx, v) with emphasis on the discrete phase-space representation of quantum tomography and 
quantum teleportation. 



III. APPLICATIONS 



Nowadays, within the context of quasiprobability distribution functions in finite-dimensional spaces, the discrete 
Wigner function has a central role in some recent researches on quantum-state tomography |l6l IT?! ] , quantum telepor- 
tation |l8l Il9t l2?l I2H. phase-space representation of q uantum computers |22j | , open quantum systems |23j ] , quantum 
information theory [2J|, and quantum computation 25]. Basically, these works are based on the well-established Woot- 
ters' approach p| for discrete Wigner functions, in which "the field of real numbers that labels the axes of continuous 
phase space is replaced by a finite field having N elements," N being the power of a prime number. Notwithstanding 
this, there are other formalisms for finite-dimensional Hilbert spaces with convenient inherent mathematical properties 
which can also be applied in the description of similar quantum systems [1, 0, IE HI El EH 13 El El El H3 ■ I n this 
section, we will show that the present formalism not only introduces new mathematical tools for the analysis of finite 
quantum systems but also can be applied, for example, to the context of quantum information processing, quantum 
tomography and quantum teleportation. 



A. Quantum information processing 



Within the most important quantum operations in quantum information processing, unitary operations have a 
prominent position [30(. Besides, in the scope of quantum information theory, the unitary depolarizers play an 
important role in quantum teleportation and quantum dense coding 0, [MJ . With respect to iV-dimensional Hilbert 
spaces, unitary depolarizers are defined on a domain Q as elements of the set 



(AO = {x e J x e x\ = xjx e = 1, e e n} 



which satisfy the relation 



i^X e OXt = Tr(0)l 



(16) 



(17) 



for any linear operator O acting on finite-dimensional vector spaces, where 1 is an identity operator. Recently, Ban 
[32| has shown that the Pegg-Barnett phase operator formalism is useful for quantum information processing as well 
as in investigating quantum optical systems. In this sense, it is worth mentioning that the symmetrized version of the 
Schwinger operator basis S(ry, £) can also be considered a unitary depolarizer, since the elements of the set 

D(JV) = {VNS(r,,0 J [V^Sfa.OHV^SfoOl 1 = [VtfS(»j,0]t[VlVSfa,O] =1, l] (18) 

obey the property 



1 

N 



iVS(77,0] T = Tr(O)l 



(19) 



This result shows that the average over all possible discrete dual coordinate and momentum shifts on the N 2 - 
dimcnsional phase space completely randomizes any quantum state defined on the finite-dimensional vector space. 
Furthermore, for s = ilj and u> <E K, the elements of the set {y/NS^'fa, (,)}%,£=-£,...,£ generalize equation HlKjl. 
being S^^rj,^) the parametrized Schwinger operator basis. Unfortunately, the implementation of such unitary 
operations in a realistic quantum-computer technology encounters an almost unsurmountable obstacle: the degra ding 
and ubiquitous decoherence due to the unavoidable coupling with the environment [33] . However, recent progress |3J] 
has developed the idea of protecting or even creating a decoherence- free subspace for processing quantum information. 
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B. Marginal distributions, Radon transforms and discrete phase-space tomography 



The marginal distributions associated with the generalized probability distribution function (fi, v) are obtained 
through the usual mathematical procedure 



z>= — £ r\——i ^ ' 



(20) 
(21) 



Note that the second equality in both definitions has been attained with the help of equation (|14|) . Consequently, the 
marginal distributions are obtained by means of discrete Fourier transforms of the s-ordered characteristic function 
calculated in specific slices of the dual plane (r), £). Now, if one considers the hierarchical process established by 
equations (|TUj) and ((TTj) . alternative expressions for the marginal distributions associated with the Wigner and Husimi 
functions can also be derived, 



q(°)(m)= E £(m'-m)Q (1 V) 



Q ( " 1) (m)= E £(m'-m)Q (0 V) 



a< >(i/)= E £(^'-^ (1 V) 

£ 

R ( ~ 1] n= E £(^-^ (0) (^') 



where the smoothing function £(x) is given by 



£(x) = 



1 i? 3 (0|ia)^ 3 (27rax|ia) + 79 4 (0|ia)i? 4 (27rax|ia) 



/2J? ^ 3 (0|ia)^ 3 (0|4ia) + i? 4 (0|ia)^ 2 (0|4ia) 



Thus, a sequential smoothing process is immediately established among the discrete marginal distributions: Q'- 1 *' (/i) — > 
Q(°'(//) — > Q( _1 )(/x) and — > !R(°)(zv) — > D^ _1 )(j/). The importance of the quantum-mechanical marginal 
distributions for s — in the context of quantum tomography in discrete phase-space has been stressed by Leonhardt 
|lfij |. where measurements on subensembles of a given quantum state are necessary in the reconstruction process. 

The Radon transforms represent an important mathematical key for quantum-state reconstruction |35j. Pursuing 
this line, Vourdas has introduced a wide class of symplectic transformations in Galois quantum systems which 
allows us to reconstruct the discrete Wigner function from measurable quantities. Basically, these symplectic trans- 
formations consist of Bogoliubov-type unitary transformations generated by J(f2i, ^2, ^3) = M(f2 3 )N(Sl2)C(ili), 
where 



c(ni) = 

N(fi 2 ) = 
M(fi 3 ) = 



1 1 

= E ex p 



N 



E ex p 

1 - 
—== y exp 



N 



(Q 2 £ - 2tj)£ 



17T 



S(r?,0) 
S(0,0 



are unitary operators written in terms of the symmetrized Schwinger basis S(rj, £), with fii = £4(1 + C2C3)~ 1 j ^2 = 
C2C4 _1 (1 + C2C3): ancl ^3 = CsC4(l + C2C3) -1 - Here, the discrete elements of the set {Cs}i=i,...,4 assume integer values in 
the closed interval [—£, £], and satisfy the relation (1(4 — C2C3 = 1 mod(iV). It is worth mentioning that this constraint 
implies in the existence of the inverse elements since (1 = C 4 T 1 (1 + C2Cs)- Now, let us initially apply the unitary 
transformation J(f2i, Q2, ^3) on the parametrized Schwinger basis 8^(77, £). Thus, after some algebra we obtain 



J(n i ,O2,^3)S (s) (^0J t ("l^2,^3) 



(22) 
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Using this auxiliary result in the calculation of J(fii, O2, il3)T( s ) (/1, i/)jt(f2i, ^2, ^3), we promptly obtain the inter- 
mediate expression 



1 

— F= / exp 



2ni 

iT&(ri,S;N)-—(w + Sv) 



s (s) (c^ + C2^C3^ + c 4 e) 



The next step consists in replacing the dummy discrete variables rj and £ by £47/ — and — ^377' in the double 
sum, respectively, with the aim of establishing the compact expression 



1 



E 



exp 



27T1 

i7r$(r?', e; N) -— (r/// + £V) 



3c(C4^-C2e,Cie-C3^) 



s«(»/,o 



(23) 



being // = — £3^ and i/ = Qiv — ^2^1 the new discrete variables written as a linear combination of the old ones. In 
particular, this transformed mod(A r ) invariant operator basis can be used to derive the marginal distributions through 
the standard mathematical procedure 



Q (s) (^;Ci,C 3 ) 
^ (S V;C2,C 4 ) 



V F {s) (u' v')8 [N] 



S [N] 



(24) 
(25) 



These results characterize the Radon transform in the present context and say that the sum of the parametrized 
function F^(/i', v') on specific lines in the ^-dimensional phase space represented by the discrete variables // and 
v 1 are equal to the marginal distributions for any value of the parameter s (when s = 0, the marginal distributions 
coincide with probabilities). In terms of the discrete s-ordered characteristic function, equations 124|) and (|25|) can be 
written as 



Q (s) (m;Ci,C 3 ) 
ft (s V;C 2 ,C 4 ) 

whose inverse expressions are given by 
S (s) (Cir?,C3?7) = 

s (s) (C 2 e,C 4 £) = 



E, 2-k'i 
ex P ( — ^-w 



■q=-t 



N 



E ex p(-?r^ 



1 

N 
1 

N 



3C(0,f) 



s (s) (Ci^,C3^) 



s (s) (C 2 ?,C 4 



DC(0,0 



E ex p 



2?ri 



]T exp(^^)^)(,;C2,C 4 ). 



w )q( s )(m;Ci,C 3 ) 



JV 



(26) 
(27) 



Note that equations i|26[l and l|27|l establish a link between measurable quantities (rhs) and discrete s-ordered char- 
acteristic functions (lhs); moreover, they can be used to construct, for instance, the quasiprobability distribution 
functions in finite-dimensional spaces. In summary, we have established a set of important theoretical results which 
constitute a discrete version of that obtained by Vogel and Risken 36] for the continuous case. 

From the theoretical point of view, the ideas of quantum computation can nowadays be used for illuminating 
some fundamental processes in quantum mechanics [37| . In this sense, Paz and co-workers |l7j have shown that 
tomography and spectroscopy are dual forms of the same quantum computation (represented by a 'scattering' circuit), 
since the state of a quantum system can be modeled on a quantum computer. Furthermore, using different versions 
of programmable gate arrays, the authors have been capable not only of evaluating the expectation value of any 
operator acting on an iV-dimensional space of states, but also of measuring other probability distribution functions 
(e.g., Husimi and Kirkwood functions) in a discrete phase-space. Here, we employ a slightly modified version of the 



9 




FT 



Measurement 




Density 
Operator 



U 



FIG. 1: Slightly modified version of the "scattering circuit" used to evaluate the real and imaginary parts of the expectation value 
Tr(Up) for a unitary operator U, where |0) represents the ancillary qubit state which acts as a probe particle in a scattering 
experiment, and H denotes a Hadamard transform. In particular, for a controlled-U operation given by U = ^/NS(r/, £), 
the measurements of the ancillary qubit polarizations along the z and y axes allow us to construct the discrete Wigner 
function W(/i, v) = TrfT 1 - ' (fi, v)p] (discrete characteristic function H*- ^ (77, £) = TrfS' ' (77, in the presence (absence) of the 
controlled-FT operation. 

scattering circuit to measure the discrete Wigner function W(/i, v). Basically, we modify this circuit by inserting a 
controlled-U operation between the Hadamard gates, with U = ^/NS(rj. £) acting on a quantum system described by 
some unknown density operator p, and also a controlled Fourier transform (FT) after the second Hadamard gate. This 
is illustrated in figure 1, where a set of measurements on the polarizations along the z and y axes of the ancillary qubit 
|0) yield the expectation values (<r z ) = v^VRe[W(/i, v)] and (a- y ) = \/^VIm[W(/i, u)], respectively. In the absence of 
the controlled-FT operation, these measurements lead us to obtain the characteristic function 2^(77, £) for s = 0, 
namely, (cr z ) = \/NRe[S^ (r], £)] and (er y ) = y/~N~kn[z,( ') (rj , £)]. However, to construct the discrete Husimi function, 
some modifications must be included in the primary circuit (see reference |l7| for more details) or the link established 
by equation between the Wigner and Husimi functions should be employed. Both situations deserve a detailed 
theoretical investigation since their operational costs can be prohibitive from the experimental point of view. Next, 
we will present a phase-space description of the process inherent to quantum teleportation for a system with an 
iV-dimensional space of states. 



In the last years, great advance has been reached in the quantum teleportation arena. In particular, we observe 
that: (i) different theoretical schemes for teleportation of quantum states involving continuous and discrete variables 
have been proposed and investigated in the literature 0, 0, |2(], El HI] , and (ii) its experimental feasibility has 
been demonstrated in simple systems through pairs of entangled photons produced by the process of parametric 
down-conversion jsfll ]. Moreover, the essential resource in both theoretical and experimental approaches is directly 
associated with the concept of entanglement, which naturally appears in quantum mechanics when the superposition 
principle is applied to composite systems. An immediate consequence of this important effect has its origin in the 
theory of quantum measurement |4fJ, since the entangled state of the multipartite system can reveal information 
about its constituent parts. 

Recently, the quasiprobability distribution functions have represented important tools in the phase-space description 
of the quantum teleportation process for a system with an TV-dimensional space of states. For instance, Koniorczyk 
et al |l9j have presented a unified approach to quantum teleportation in arbitrary dimensions based on the Wigner- 
function formalism, where the finite- and infinite-dimensional cases can be treated in a conceptually uniform way. 
Paz |2(j has extended the results obtained by Koniorczyk et al to the case where the space of states has arbitrary 
dimensionality. To this end, the author has used a different definition for the discrete Wigner function which permits us 
to analyze situations where entanglement among subsystems of arbitrary dimensionality is an important issue. Here, 
we use the new mod(iV ) invariant operator basis TW (^. 1/) i n order to obtain a discrete phase-space representation of 
quantum teleportation which permits us to extend the results reached by Paz in the discrete Wigner-function context 
for any discrete quasiprobability distribution functions. 



C. Discrete phase-space representation of quantum teleportation 
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Generalized Bell states 



The generalized Bell states were first introduced by Bennett et al |18| in the study of quantum teleportation for 



systems with N > 2 orthogonal states. Basically, these states can be defined as , 



Vf^TV 2 |*o,o), where 



1 1 

l*o,o) = -j= l«e)l®K)2 



represents the pure state maximally entangled for a bipartite system (in this case, the reduced density matrix of each 
constituent part is equal to (1/A r )l i for % = 1,2), being {\v a )i} a = -£,..., I and {|tt^)j}/9=-4,...,^ the eigenstates of the 
Schwinger unitary operators Vj and Uj, respectively. Furthermore, the generalized Bell states satisfy the following 
properties: 

(i) (^uAjm) =^/l> 2 (orthogonality relation) 

t 

(ii) ^2 \^u>i,u 2 ) (*wi,«2 1 = !i ® I2 (identity relation) 

(iii) U+l*^,^) =U 1 ®U 2 |*c 1 , W2 ) = exp[-(2 7 ri/JV0«i]|¥ Wl 
V_|* Wl>£l)2 ) = ViigiVa 1 !*^) =ex P [(27ri/JV)w 2 ]|* Wl , 



W 2 ) 



Note that U+ displaces both systems in coordinate by the same amount u)\, while V_ displaces them in momentum 
by the quantity u> 2 in the opposite direction. In addition, as {|vE , WljW2 )} Wl!W2= -£,... l £ are common eigenstates of U + and 
V_ , such states can be interpreted as corresponding to the eigenstates of the total momentum and relative coordinate 
operators I 19|. l20l : indeed, these states are the discrete version of the continuous ones used by Einstein, Podolsky, 
and Rosen [4l|. Thus, the generalized Bell measurements will be characterized in our context by the set of diagonal 

projection operators {|* Wl)Ua )(*u<i,w2l}wi,w a =-i,...,«- 

Now, let us establish some further results related to the generalized Bell states and their discrete phase-space repre- 
sentation. The first one corresponds to the mapping of \^>ui.u 2 ) I m terms of the basis {T^ s, ^(/Uj, Vi)} ^ uUi =-(.,...,i 
for each subsystem, i.e., 



1 1 

I^ 2 )<*^ 2 | = w rt-' 1 '-"Hui,<*MMW,*,*i,^)Tr ) {Mi,'*)' 



N 2 

fJ.l,Vl,fl2,V2 =—i 

with the coefficients of the expansion given by 



(28) 



Consequently, the second one refers to the inverse mapping of l)28[l. which can be directly reached with the help of 
property (ii) as follows: 



T^Vi^i)®Tt 2 V 2 ^ 2 )= Y e^^)^!,^,^,^!^,^,^,^)!^,^)^^! (29) 

a;i,aj 2 ,a;^,aj 2 = — i 



being 



Q^> s *\ f i 1 ,v 1 ,Li 2 ,v 2 \w 1 ,u> 2 ,<j' 1 ,u>' 2 ) = Tr T^Vi^i) ® T^ a) (^ 2 , ^I^.^X* Wl 



It is worth mentioning that a general connection between the coefficients of both expansions (|28|) and l|29|) can also 
be promptly established for any values of {s\, s 2 } € K, 

T (_Sl '^ S2 '>(u;i,LJ 2 ,Uj' 1 ,LJ 2 \fJ,i,V 1 ,fX2,V2) = e (_Sl ' _S2) (^l,^l,Ai2,^2ki,W 2 ,Wi,W 2 ) 

The analytical expression of these coefficients will be omitted here due to its apparent irrelevance in the phase-space 
description of the quantum teleportation process. 
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However, some useful results derived from these coefficients deserve to be mentioned and discussed in detail. For 
instance, equation l|29[) allows us to calculate the parametrized function 

J&'V.^/tt.^Hn^ (30) 

which coincides with for particular values of u>i and u>[. In this situation, the analytical expression 

^ 2 2) (Mi,n,/i2^ 2 ) = ^ E ex P [^[^ 1+M2 + Wl ) + ^i-^-^)]}[^,er (si+S2) 

can be reduced to the following discrete quasiprobability distribution functions: 

• si = s 2 = (Wigner function) W wl)W2 Oi,!/i,/x 2 , fa) = ^_( Atl+ ^ 2 )^ 1 / 1 -^ 



si — s 2 — -1 (Husimi function) ^ 1 ,ui 2 (^i,i'i, ^2,^2) = N 1 \3C((ii + /i2 + wi, ^1 - v 2 - w 2 )| 2 



To measure the discrete Wigner function associated with the generalized Bell states, some minor modifications should 
be implemented in the scattering circuit (see figure 1): the first one concerns to the controlled-U operation between the 
Hadamard gates, since it must be replaced by U = [yiVSi^i, £1)] <g> [VNS2{t)2, £2)} in order to process operations for 
bipartite systems; while the second one consists in preparing the input density operator in the generalized Bell states, 
namely p = l^wi.waK^wi.u^l' This procedure leads us to obtain the expectation value (<r z ) — NW Ul;U12 (fjti, fx, (12, ^2) 
through a set of measurements on the polarization along the z-axis of the ancillary qubit. Furthermore, these minor 
modifications on the scattering circuit can also be used to measure any discrete Wigner function associated with a 
general bipartite system. 



2. Quantum teleportation 



Basically, the quantum teleportation process consists in a sequence of events that allows us to transfer the quan tum 
state of a particle onto another particle through an essential feature of quantum mechanics: entanglement |l8t l39|. 
In this sense, let us introduce a tripartite system described by p = p 1 ® (|^E f o,o)(^ / o.o|)23, where subsystems 2 and 3 
were initially prepared in one of the Bell states. The plan is to teleport the initial state of subsystem 1 through the 
protocol established in [20| . 

1. We initiate the protocol considering the density operator associated with the tripartite system written in terms 



of the new basis {T\ Si '(fii, 



for each subsystem i = 1, 2, 3 as follows: 



P = 

where 



1 

N3 



E 



( — 52, — S3) 



(M2, V2, A*a, ^)T< si) (mi, vi) ® T% 2 \n2, ^2) <8> T^Qa, v 3 ) 



Ml ,^1,^2,^2,^3, V 3 = — l 



F[ Sl \^ Vl ) = Tn \t[ S1 \ i x 1 ,v l )p 1 



r-,(-s 2 ,-s 3 ) 
23 



(M2, ^2,^3,^3) = Tr 



2:5 



(^2,^2) 



,(-53) 



(At 3 ,^)(|*o,o)(*o,o|) 



23 



Next, we perform a measurement on subsystems 1 and 2 that projects them into the Bell states (this procedure 
corresponds to a collective measurement which determines the total momentum and relative coordinate for 
composite subsystem 1-2). For convenience, before the generalized Bell measurement, let us express the phase- 
space operators (/ii, v{) (g) T 2 {^2, V2) according to equation 1291) . 



1 : 

7V3 



Ml, l/ l,M2,^2,A 1 3,^3 = -^ 

E e (si ' S2) (^i,^i,^ 2 ,^2|wi,W 2 ,wi,W 2 ) (l*a) ljWa }(*u4 jtl > 2 |)l2 <H> T^ 3 ] (/i 3 , I*) 
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Thus, after the measurement on the first two subsystems, only the terms with lo\ = u>[ = a and u>2 — lo' 2 — /3 
survive. Consequently, a reduced density operator for the third subsystem can be promptly obtained, 



PSR=^ E A^ 1 - S3 )(m3^ 3 )T| S3) (M3^3), (31) 

which does not depend on the complex parameter S2- Here, the coefficients are given by 



with 



<r i) (^^i^3^ 3 ) = ^ E cxpl^K^-Ma + a)-^!-^-^)]}^^)] 33 ^ 1 • 

Note that p 1 ' S3 \fj,^, 1/3) simply tells us how to construct, independently of parameters si and S3, the final 
parametrized function for the third subsystem from the initial parametrized function of the first one. 

2. Now, let us analyze the particular case si — S3 = s. In this situation, equation (|31|l assumes the simplified form 

1 1 

P™=N E ^i ( ~ S) (M3-a^3+/3)TfV 3 ^3) (32) 

where the coefficients F± ^ ([13 — a, v% + /3) play a central role in the phase-space description of the quantum 
teleportation process. In fact, they allow us to conclude that, after the generalized Bell measurement, the third 
subsystem has a parametrized function which is displaced in phase-space by an amount (— a, (3) with respect to 

the initial state of the first subsystem, namely F^ R S \^, v^) = f[ (/is — a, v% + /?). Therefore, the recovery 
operation basically depends on the calibration process of the generalized Bell measurements performed on the 
first two subsystems: for instance, when a — (3 — 0, we reach a complete recovery operation. 

In short, we have presented a quantum teleportation protocol that leads us to obtain a phase-space description 
of this process for any discrete quasiprobability distribution functions associated with physical systems described by 
iV-dimensional space of states. 



IV. CONCLUSIONS 



In this paper we have employed the new mod(iV) invariant operator basis {T( S J (/i, v)} n,v=-i,...,i recently proposed in 
|26| , with the aim of obtaining s-parametrized phase-space functions which are responsible for the mapping of bounded 
operators, acting on a finite-dimensional Hilbert space, on their discrete representatives in an iV 2 -dimensional phase 
space. In fact, we have established a set of important formal results that allows us to reach a discrete analog of the 
continuous one developed by Cahill and Glauber [jj. As a consequence, the discrete Glauber-Sudarshan (s = 1), 
Wigner (s = 0), and Husimi (s = —1) functions emerge from this formalism as specific cases of s-parametrized 
phase-space functions describing density operators associated with physical systems whose space of states has a 
finite dimension. In addition, we have also established a hierarchical order among them that consists of a well- 
defined smoothing process where, in particular, the kernel 3C(?y,£) performs a central role. Next, we have applied our 
formalism to the context of quantum information processing, quantum tomography, and quantum teleportation in 
order to obtain a phase-space description of some topics related to unitary depolarizers, discrete Radon transforms, 
and generalized Bell states. Indeed such descriptions have allowed us to attain new important results, within which 
some deserve to be mentioned: (i) we have shown that the symmetrized version of the Schwinger operator basis 
{S(ry, C)}j7,f=-^,...,f can be considered a unitary depolarizer; (ii) we have also established a link between measurable 
quantities and discrete s-ordered characteristic functions with the help of Radon transforms, which can be used to 
construct any quasiprobability distribution functions in finite-dimensional spaces; and finally, (iii) we have presented a 
quantum teleportation protocol that leads us to obtain a generalized phase-space description of this important process 
in physics. It is worth mentioning that the mathematical formalism developed here opens new possibilities of future 
investigations in similar physical systems or m the study of dissipative systems, where the decoherence effect has 
a central role in the quantum information processing (e.g., see reference 23]). These considerations are under current 
research and will be published elsewhere. 
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